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Abstract 

A microscopic description of the zero energy two-body ground state and many-body static prop- 
erties of anisotropic homogeneous gases of bosonic dipoles in two dimensions at low densities is 
presented and discussed. By changing the polarization angle with respect to the plane, we study 
the impact of the anisotropy, present in the dipole-dipole interaction, on the energy per particle, 
comparing the results with mean field predictions. We restrict the analysis to the regime where the 
interaction is always repulsive, although the strength of the repulsion depends on the orientation 
with respect to the polarization field. We present a series expansion of the solution of the zero 
energy two-body problem which allows us to find the scattering length of the interaction and to 
build a suitable Jastrow factor that we use as a trial wave function for both a variational and 
diffusion Monte Carlo simulation of the infinite system. We find that the anisotropy has an almost 
negligible impact on the ground state properties of the many-body system in the universal regime 
where the scattering length governs the physics of the system. We also show that scaling in the 
gas parameter persists in the dipolar case up to values where other isotropic interactions with the 
same scattering length yield different predictions. 
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I. INTRODUCTION 



Quantum dipolar systems of bosons and fermions have gathered much experimental and 
theoretical attention in recent years. In 2005, Griesmaier and collaborators [l] on one 



side, and Stuhler and collaborators p| on the other, reported on the first experimental 
realization of a Bose condensate of ^^Cr, where the dipolar moment of the atoms is so large 
(~ Qfis) that the effect of the dipole-dipole interaction is comparable in strength to the van 
der Waals forces. More recently, new and exciting results have been achieved with polar 
molecules of Rubidium and Potassium f'^^K^^Rb) |3|, which have not been easy to create 
due to strong loss rates in the population induced by chemical reactions [5|, |6|. A promising 
route towards a molecular Bose-Einstein condensate is Feshbach association of Rb and Cs, 
which are not reactive j4|. One of the major advantages of polar molecules is that the electric 
dipole moments are remarkably larger than in the magnetic case of ^^Cr and can be tuned 
by applying an external electric field. Systems of polar molecules have been speculated to 
present interesting applications ranging from the control of chemical reactions 7^ to practical 



applications of quantum information theory 
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From the theoretical point of view, dipolar systems present novel and interesting phenom- 
ena that make them particularly appealing. On the one hand, the anisotropic character of 
the dipole-dipole interaction introduces additional degrees of freedom compared with other 
condensed matter systems that can potentially enrich the phase diagram. On the other, the 
interaction decreases at large distances as and becomes long ranged in three dimensions 
(3D), in contrast to typical van der Waals forces. In two dimensions (2D), though, the 
interaction is still short ranged but at the border between both regimes. 

The potential Vd{r) describing how two dipoles with dipolar moments pi and p2 interact 
is given by 

(1) 



Pi ■ P2 - 3(pi ■ r)(p2 ■ f) 



with r the relative position vector between them and Cdd the coupling constant defining the 
strength of the interaction. For permanent magnetic dipoles Cdd = /Wo/^^ where /iq is the 
permeability of vacuum and fi is the permanent dipole moment of the atoms. Alternatively, 
the electric dipole moment can be induced by an electric field E, and in this case the coupling 
constant is Cdd = d'^/^o, where d = aE with a the static polarizability and eo the permitivity 
of vacuum. For a system of fully polarized dipoles in 2D as the ones considered here, pi and 
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P2 are parallel and define a fixed direction in space, see Fig. [T] In this case Vd{r) simplifies 
to 



VJr) 



a 



dd 



1 - 3A2 cos2 e 



(2) 



where A = sin a, a being the angle formed by the normal to the plane and the polarization 
field, which is tilted towards the x-axis. In this expression, r and 6 stand for the in-plane 
distance and polar angle, respectively. Notice that, in contrast to what happens in three 
dimensions, a is fixed in the fully polarized system and thus A < 1 is a constant of the 
problem for a given a. 

One of the consequences of the form of the interaction is that it can be either attractive or 
repulsive depending on the orientation with respect to the polarization field. This fact alone 
triggers interesting discussions about the static and dynamic properties of dipolar systems. 
For instance, the formation of a density instability observed in the elementary excitation 
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spectrum as the formation of a deep roton minimum has been widely discussed 
feature that is absent when the interaction is reduced to its purely isotropic limit in 2D 
corresponding to A = in Eq. The presence of two- and/or many-body bound states 
in stacks of dipolar layers has also raised interesting questions [13]. Many other aspects 



including scattering properties in 2D |14j4l6|. pseudopotential treatments [l7| or the impact 



of the anisotropy on the superfiuid properties of a dipolar Bose gas [18^ have gathered interest 
in the recent years. Most of these aspects and many others are covered in detail in the review 
article by Lahaye et al 19|] . In this work we concentrate on the case where the interaction is 
always repulsive, but is still anisotropic. That defines a critical angle ttc = 0.615 above which 




FIG. 1: Two dipoles confined to move on the X-Y plane. The polarization field lays on the XZ 
plane and fixes a direction in space forming an angle a with the z axis. 
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Vd{r) starts to present negative contributions. We thus analyze the role of the anisotropy of 
the interaction in situations where the system is clearly prevented from collapse. Since etc 
is reached at = 1/3, we propose to use as an expansion parameter. 

Despite the relevance of the anisotropic character of the interaction in all these effects, 
not much effort has been put in the microscopic description of the wave function. The 
anisotropy makes the different angular momentum channels couple in a non-trivial way, and 
that introduces additional degrees of freedom that must be taken into account in a numerical 
simulation. At low densities, though, the problem simplifies since the zero-energy solution 
of the two-body scattering problem dominates the ground state many-body wave function 
of the homogeneous gaseous phase. In this work we provide a detailed description of the 
zero-energy scattering wave function corresponding to two dipoles moving on the plane, as 
a function of the polarization angle. Equipped with that solution we build a microscopic 
variational many-body wave function that we use in a Monte Carlo simulation to describe 
the static properties of a gas of polarized dipoles in 2D at low densities. 

The analysis of the low density equation of state of a gas of weakly interacting particles 
has historically attracted great interest. Corrections to the mean field prediction for three- 



dimensional I20il and one-dimensional 



2l| systems have been known for more that fifty years 



now. The two-dimensional case has been much more controversial as already the two-body 
problem presents logarithmic divergences in the leading scattering parameters that make 
series expansions difficult to carry out 22|, |23|. In any case, the low density behavior of 
a gas of weakly interacting particles in 2D has been widely discussed in the literature for 
the case of isotropic interactions. One of the most remarkable properties exhibited by these 
systems is the universal behavior of the energy per particle, which admits a non-analytic 
series expansion in the gas parameter x = no?, with n the density and a the s— wave 
scattering length. The leading order, mean-field term in this series has been derived by 



several authors 
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25[ and reads 



( 



2'mo? 



E 

N 



Attx 



(3) 
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The detailed form of the next-to-leading correction to this expression have been the subject 
of discussion and different authors proposed different forms in the past, see for instance 



Refs. 
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271]. The correct expression was recently derived in Ref. 



numerically intense Monte Carlo simulations in Ref. 



and checked against 



29| . For the model system of hard 



disks, the mean field prediction of Eq. ([3]) fiolds well starting at x ~ 0.001 and down to 
quite low but still experimentally affordable values of the gas parameter 0]. However, 
no particular attention has been paid in all these works to the special case of anisotropic 
interactions. 

In this article we discuss to which extent the mean field law of Eq. holds for the special 
case of the spatially anisotropic dipole-dipole interaction of Eq. ([2]) when the polarization 
angle varies between and ac- We are particularly interested in discerning whether the 
angular dependence of the interaction has a noticeable impact on the mean-field prediction 
of Eq. ([3]) and on other relevant ground state properties. In order to do that, we first solve 
in section II the zero-energy two-body scattering problem and obtain an expression for the 
scattering length as a function of the polarization angle. We then use this result to build 
in section III a variational many-body wave function of the Jastrow form that we use as 
an input to both a variational Monte Carlo (VMC) and a diffusion Monte Carlo (DMC) 
calculations from where we obtain the equation of state as a function of x. For the sake of 
completeness we also analyze the pair distribution function, the static structure factor, the 
one-body density matrix and the condensate fraction, and discuss how these quantities scale 
on the gas parameter for different polarization angles and densities. Finally, in section IV 
the main conclusions of the work are summarized and discussed. 
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II. ZERO ENERGY TWO BODY PROBLEM 



In this section we develop a series expansion of the zero energy scattering solution for 
two dipoles moving on the XY plane as a function of the polarization angle a. This is done 
by first building the Green's function of the a = case corresponding to a polarization field 
perpendicular to the plane. This is a particularly suitable situation since in that case the 
interaction is isotropic (A = in Eq. ([2])) and the Schrodinger equation can be exactly solved. 
From there we carry out a series expansion in powers of of the solution corresponding 
to the anisotropic case. The anisotropy makes the different angular momentum channels 
couple, and we report the expression for all orders in A^ contributing to each partial wave. 
We end this section reporting the scattering length of the dipole-dipole interaction as a 
function of the polarization angle, required afterwards to analyze the low density properties 
of the anisotropic many-body system. 

The Hamiltonian describing the relative motion of two polarized dipoles of mass m moving 
on the plane reads 



2M 47r 



l-3A"cos2^1 

(4) 



where M = m/2 is the reduced mass. In the following we use dimensionless variables scaled 
according to the characteristic dipolar length = mCdd/ ^t^^ and energy = /mr'\. 

The two linearly independent solutions of the zero energy Schrodinger equation (SE) for 
the isotropic case are the building blocks we need in order to generate the Green's function 
that we will use afterwards to solve the anisotropic case. Setting A = 0, the SE of the 
relative motion of the two dipoles at zero energy reduces to 
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^ + ^¥^ = 0, (5) 



and the general solution of this equation can be expanded in partial waves as 

oo 

'fir, 0) = ^ ifnir) cos(n^) , (6) 



n=0 



where the wave function of each separate mode fni'f') satisfies 



''^^%Vfi + 3l»'»-0. (7) 



r dr \ dr J \ ^ 
which is a modified Bessel equation for ip„(2/\/r). The two linearly independent solutions 



of Eq. ([7]) are K2n{'2/y/r) and /2n(2/\/r) |3l|, and these two functions enter in the Green's 



function we write below. The zero energy solution of the SE requires the condition ip2n{f = 0) 
to be zero, and that discards the /2n(2/-\/r) contributions. Apart from a normalization 
constant, the zero energy solution of the A = problem becomes then 

^nir) = K,J^] . (8) 



The physical solution of the isotropic case corresponds to n = as otherwise anisotropic 
contributions would dominate at large distances since fC2n(2/-\A') grows as r" when r — > oo. 
In ID and 3D, the knowledge of the E = solution allows one to immediately obtain an 
exact expression for the s— wave scattering length a^. The description of two-body scattering 
in 2D is more involved since the low-energy expansion of the scattering amplitude diverges at 
low energies, thus introducing additional problems not found in higher and lower dimensions. 
The scattering length can however be defined to be equal to the position of the node of 
the asymptotic form of the zero energy two-body wave function. This definition has the 
additional advantage that can be used in any dimensions and will therefore be adopted 



throughout this work 
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32l |. The large r behavior of Ko{2/y/r) is —7 + | ln(r) and that 



yields the well known expression 

a, = 6^^^ 3.17222... (9) 

where 7 is Euler's gamma constant. 

The SE describing the anisotropic case can be cast in the form 

-vV + ;:^0 = ^:^0, (10) 

and the general solution for A 7^ can be derived from the Green's function corresponding 
to A = 0, which fulfills the equation 



+ G(r,rO=<5(r-r') , 



(11) 



leading to 

cos^ ft 



r cos"' f) 

0(r) = Mr) + SA^ J dy —^G{r, y)0(y) (12) 

with (foir) = Ko{2/^/r) the A = solution as described above. 

Equation ( ITTl) can be solved expanding the Green's function in partial waves as before 

1 1 °° 

Gir, r') = —goir, r') + - ^ fi'„(r, r') cos [n{e - 9')] (13) 



27r TT 

n=l 



where 

9nM=\ f (14) 

I (^) (^) if r > r' 
satisfies the boundary condition gn{r = 0,r') = gn{r,r' = 0) = while keeping it bounded 
at large distances. The general solution of the Fredholm integral equation ( 1121) admits a 
series expansion in powers of 



oo 



v^vM-z^AV^Hr), (15) 

k=0 

where each (p^'^^r) satisfies the recurrence relation 

<P^'^'\r) = 3X'I rfy ^^G(r,y)#)(y) . (16) 

When </)'^^)(r) is further expanded in partial waves and the Bose symmetry is taken into 
account 

oo 

<P^'\r,e) = Y,4nir)cos{2ne) , (17) 

n=0 

the coupling between the different angular momentum channels produced by the cos^ 6 term 
of the interaction emerges and the radial functions satisfy the following recurrence relations 
for even n 

€^'\r) = ^-^fdy ^ i^nUy) + 2<pL'\y) + <p\tUy^] (is) 

that can be solved iteratively starting from 0o°^(r) = (po{r). From this expression one 
sees that by adding successive orders in to the series expansion of 0(r), more angular 
momentum channels couple together. As in the regime considered the interaction is fully 
repulsive, A < l/y/S and that makes A^ a small parameter that we can use in a series 
expansion of the solution. In fact, it can be shown from the previous expressions that 
(f)^n{r) = for 2n > k, and that therefore the lowest order contribution to the n-th mode is 
A2". By adding 

02n ) for all k and fixed n one recovers 02n i^" ) ; the complete 2n — th mode 
contribution to 4>{r). We thus find 



oo 

02n(r) = J^A^VSW 

k=n 



which means that, up to a given order A^^, the total wave function 0(r) has contributions 
coming only from channels n = 0,2, ... ,2k. 
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FIG. 2: (Color online) Lowest order partial wave contributions to the zero energy two-body scat- 
tering wave function. The upper and lower panels show (j)^^^ = Kq{2/^) (black dotted line), (j)o{r) 
(red solid line), 02 (^) (green dashed line) and 04 (r) (blue dot-dashed line) for the polarization 
angles a = 0.2 and a = 0.6. The inset show the a = 0.6 curves in logarithmic scale. 

Figure [2] shows the lowest order partial waves contributions corresponding to the two 
polarization angles a = 0.2 (upper panel) and a = 0.6 (lower panel), the latter being the 
largest angle considered in this work. The latter angle is fairly close to the critical angle 
Oc = 0.615 where the interaction ceases to be purely repulsive. At a = 0.6 the contribution 
of the n ^ modes is expected to be larger than for any lower angle. This is clearly seen 
from the figure, where the corrections to the n = and n = 2 partial waves are shown 
(red solid and green dashed lines), as well as the leading correction corresponding to 
the n = mode. It is clear from Eqs. ( IT4l) to ( fT8|) and the positiveness of the modified 
Bessel functions that every radial contribution (j)^2n l'") ^'^ two-body wave function is also 
positive, as seen for the lowest mode contributions in the figure. It is also apparent that the 
lower the angle, the smaller the correction to the a = solution (p{r) is, as expected. Despite 
the fact that the series expansion of the two-body solution 0(r) is in general alternating due 
to the cosine terms, the total two-body wave function does not change sign as the interaction 
is everywhere repulsive, thus making the E = scattering solution be the ground state. 




9 




FIG. 3: (Color online) Cuts of the zero energy two-body scattering wave function describing the 
relative motion of two dipoles. The blue solid line and the red dashed line correspond to the cuts 
along the x and y axes, respectively. The green circles and black squares are the prediction of the 
optimal Jastrow factor obtained from the solution of the HNC / Euler equations for a value of the 
gas parameter x = 0.01. 

The effect of the anisotropy on the ground state wave function is seen in Fig. [3] where 
two cuts, one along the x-axis (contained in the plane formed by the polarization vector 
and the 2;-axis), and another in the perpendicular direction (y-axis) are depicted for the two 
polarization angles a = 0.2 and a = 0.6. These cuts coincide with the directions where the 
interaction is least and most repulsive, respectively. As it can be seen, anisotropic effects 
are visible in both cases but are more pronounced at high polarization angles . The inset in 
Fig. [2] shows that the m = mode dominates at very large distances as expected, making 
the asymptotic wave function be isotropic. From (f)Q{r) one can extract the scattering length 
as{X) of the anisotropic dipolar interaction, which is given by the node of its asymptotic 
r — )■ oo form. An analytic approximation to as(A) can be easily obtained recalling that every 
mode 0„(r) contributes to order A^" and that therefore the anisotropy enters at order A^. 
Direct inspection of the modes expansion of the the Schrodinger equation for (f){r) reveals 
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that (poi'f') and (f)2{f) are related according to 



1 d ( rf0o\ 1 , 3A 
1 + —00 



2 r 

2r3 



^^^r) + -02(r) 



(19) 



r dr \ dr J 

An approximation of order to as(A) can be obtained by keeping only the (poi^) niode on 
the right hand side of this equation 




1 d 
r dr 

which once again is a modified Bessel equation with the general solution 



(20) 





3A£ 
2 



(21) 



with Af a normalization constant. In this expression -B(A) is an unknown function of A^ 
according to the parity of the Hamiltonian under the A — )■ — A transformation. Furthermore, 
-B(O) = so that one recovers the isotropic solution given in Eq. (Q. Hence -B(A) = b2X^ to 
order A^, with 62 a constant. In the asymptotic r — )■ 00 regime, /q ~ 1 and one can compare 
the expansion to order A^ of the above expression to the expression of (f)o{r) to the same 
order obtained from the integration of the Green's function done before. This yields 62 = 
and one has 

1 



Mr ^00)^ Ko 1^21 
with as(A) the s-wave scattering length 

a«(A) = 



3Af 

2 



In 7TT 

2 a,(A) 



(22) 



,27 



3A^ 
2 



(23) 



This expression is accurate up to order A^, so one could expect it to provide a reliable 
prediction only at small polarization angles. This turns out not to be the case, and in Fig. H] 
we show the comparison of this approximation to the exact result obtained by numerically 
finding the node of the asymptotic m = wave function, which is isotropic and dominates 
the large distance behavior of 0(r). As can be seen from the figure, the approximation works 
surprisingly well up to the critical angle etc where the interaction ceases to be fully repulsive. 
Deviations increase with increasing polarization angle, but even at a = the separation 
between the approximation in Eq. ( 123|) and the exact numerical estimation is less than a 
3%. 
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FIG. 4: (Color online) s-wave scattering length as(A = sin(Q)) as a function of the polarization 
angle a. The blue solid line and the black dots correspond to the exact numerical solution and the 
approximation of Eq. ([23|) . respectively. 

III. MANY-BODY DESCRIPTION 

In this section we extend the previous discussion and analyze, using DMC and VMC 
methods, the most relevant ground state properties of an homogeneous and anisotropic gas of 
bosonic dipoles with a polarization angle a lower than the critical angle ac = 0.615. We stick 
to the low density limit were the system, characterized by a fully repulsive and anisotropic 
interaction, remains in a stable, gaseous phase. One of the most relevant quantities to 
analyze at low densities is the total energy per particle of the gas and its universal scaling 
properties. Quite a lot of work has been devoted in the past to that question, including both 



3D 



20 



33 



35| and 2D |24j-l30l| systems. However, little has been discussed about the same 



properties in anisotropic systems as the dipolar gas considered here. We analyze the impact 
of the polarization angle a formed by the dipoles on the universality scaling law exhibited 
by other isotropic, short ranged interactions. 

The Hamiltonian of the system of fully polarized dipoles, written in the same dipolar 
units used in the previous section, becomes 



1 - 3A2 cos2 a 



(24) 



j=i i<j u 

with A = sin a, and fj, and 9ij the distance and angle formed by dipoles i and j, respectively, 
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measured on the plane. 

The leading ground state quantities describing the low density static properties of the 
system can be obtained using different techniques. In this work we stick to diffusion and vari- 
ational Monte Carlo methods, widely used nowadays in the analysis of weakly and strongly 
correlated systems. Variational Monte Carlo samples stochastically a trial wave function 
and yields an upper bound to the real ground state energy of the system. On the other 
hand, diffusion Monte Carlo uses also a trial wave function to guide the sampling process 
but removes the contributions from excited states to finally yield statistically exact energies. 
In both cases, however, a suitable variational wave function is required. The quality of the 
results is directly related to the quality of the wave function employed in the VMC case, 
while DMC is far less demanding and any reasonable guiding function can be used as long 
as it is not orthogonal to the true ground state. But even in DMC a high quality wave 
function makes the method converge faster and with smaller variance towards the exact 
result. Consequently, seeking for a good trial many-body wave function \l/(ri, r2, . . . , r^) is 
always desirable. 

In this work we use a model wave function of the Jastrow form 

vl/(ri,r2,...,r;v) = n/(r^i) ' (25) 

i<j 

where the two-body correlation factor f{rij) = f{ri — rj) depends on the position vector 
linking particles i and j. One significant difference between this Jastrow factor and the ones 
most commonly employed in the analysis of other condensed matter systems is that, due to 
the anisotropic character of the interaction, /(r) depends explicitly on the whole r vector 
and not only on its magnitude. In this way, the wave function in Eq. (I25ll describes an 
homogeneous but anisotropic system as the one under study. 

At low densities, the zero-energy scattering solution of the two-body problem greatly in- 
fluences the structural properties of the gas. For that reason we use as a Jastrow factor the 
anisotropic solution of the relative motion of two dipoles on the plane derived in the previous 
section. The n = mode of this wave function is matched at some healing distance ^ with 
the symmetrized form of a phononic wave function /^(r) = exp(— C/r) 36j, taking both ^ 
and C as variational parameters and imposing the continuity of /(r) and f'{r) at r = ^. 
The n > modes of the two-body problem, inducing the anisotropy of \l/(ri, r2, . . . , r^r) in 
Eq. (|25|1 . decay to zero at large distances and so their influence at the boundaries of the 
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TABLE I: DMC and VMC energies per particle as a function of the gas parameter x = nc? . 



simulation box is marginally small. Alternatively, the optimal Jastrow factor corresponding 
to the many-body problem can be obtained from the solution of the HNC/0 Euler-Lagrange 



equations 



37|. Although not exact, the optimized HNC/0 solution gives an accurate varia- 



tional description of quantum Bose systems and captures most of the short and long range 
features of the exact ground state wave function. For the sake of comparison, we also show 
in Fig. |3]the optimized HNC/0 Jastrow factor (black and green symbols) at a; = 0.01 and 
polarization angle a = 0.6. The comparison indicates that the two-body solution provides 
an accurate description of the two-body correlation factor, which becomes even better as the 
gas parameter is reduced. We have checked that the HNC/0 Jastrow factor and the solution 
of the two-body problem are in very good agreement in the whole range of gas parameter 
values considered in this work. 

Table H] lists both the VMC and DMC energies obtained from the Jastrow trial wave 
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FIG. 5: (Color online) Ratio of the energy per particle of the gas of dipoles for different polarization 
angles to the mean field prediction of Eq. ([3]). Black circles, red triangles and blue triangles 
correspond to q = 0.2,0.4 and 0^, respectively. The green squares are the optimized HNC/EL 



energies for hard disks of Ref. 



35| |. while the solid line is the universal curve of Ref. 



33]. The 



dotted line corresponds to the mean field prediction. 

function of Eq. f l25p for the polarization angles a = 0.2, 0.4 and 0.6. Notice that the energies 
in the table are given for fixed x and different polarization angles, and since the scattering 
length varies with a, the densities change accordingly. A direct measure of the quality of the 
variational model is given by the separation between these two measures (VMC and DMC), 
and one can check that the relative difference in energies is always of the order of 1% or 
2%. Other than that, the energy is an increasing function of the gas parameter that yields 
appreciably different results for different polarization angles. These energies can be used to 
check the influence of the anisotropic character of the dipolar interaction on the universality 
scaling property fulfilled by the energy per particle of homogeneous and isotropic systems 
in 2D. In order to do that, one has to express the total energy per particle in units of 
li^/2mo? with a the scattering length. This is achieved multiplying the energies in Table H] 
(expressed in dipolar units) by 2a^(A), with as{X) the scattering length for the corresponding 
polarization angle. Figure shows the ratio of the energy per particle in units of f? j2mc? 
to the mean field prediction of Eq. (|3]) for the three polarization angles a = 0.0, 0.4 and 0.6. 
As it can be seen, expressed in scattering length units, all curves corresponding to different 
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polarization angles merge into a single curve, with very small deviations that are not easily 
resolved even at the highest values of gas parameters x considered in this work. That means 
that the anisotropy of the interaction, present in the wave function, does not appreciably 
affect the energy per particle in the low density regime analyzed in this work. We conclude 
that the difference in energy values shown in Table [I] for fixed x and varying polarization 
angles are to be mostly attributed to the different density n = x/a? in each case. 

Figure [5] also shows the universal curve including beyond mean field effects of Ref. 29 1 
and the optimized HNC/0 prediction for a gas of hard disks of Ref. 30[]. As it can be seen, 
the universal and the hard disks curves are close to each other while the dipole curves remain 
closer to the mean field prediction 6^/(2;) as the gas parameter is raised. Starting at x ~ 0.05 
the dipole curves bend downwards and the energy deviates significantly from emf{x). In any 
case, it is clear from the figure that the universality regime where the energy per particle 
depends only on the gas parameter of the interaction is left much before anisotropic effects 
have an appreciable impact on the energy of the dipolar gas. 

The anisotropic character of the dipolar interaction has a direct influence on the ground 
state wave function that is reflected in the ground state expectation value of any many-body 
operator. Figures |6] and [7] show pure DMC estimations [38] of the pair distribution func- 
tion (yf(r) and its Fourier transform, the static structure factor S'(k), for two values of the 
polarization angle a = 0.4 and a = 0.6 (left and right panels), and three values of the gas 
parameter x = 10^^, 10^^ and 10^^ (top to bottom). Notice that in both figures the horizon- 
tal axis has been scaled with the square root of the density for a better comparison. Due to 
the symmetries of the Hamiltonian, the complete ^'(r) and >S'(k) functions vary continuously 
on the plane but the pattern on the first quadrant is repeated and reflected on the other 
three. The figures show only the two cuts along the perpendicular and parallel directions 
with respect to the polarization plane, corresponding to the lines where the interaction is 
most and least repulsive, respectively. As it can be seen, and in agreement with what one 
would expect, the effect of the anisotropy is more clearly seen at higher polarization angles 
and for large values of the gas parameter, being maximal for a = 0.6 and x = 10^^. For 
fixed a the separation between g{r, 0) and g{0, r) is enhanced with increasing x, as happens 
with S{k, 0) and S{0, k). Accordingly and for a given x, the separation between the curves 
also increases when the polarization angle is raised. In any case it is remarkable how the 
anisotropy present in g{r) and 5'(k) changes with the polarization angle as can be seen from 
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FIG. 6: (Color online) Pair distribution function for a = 0.4 and a = 0.6 and three values of the 
gas parameter. The red and blue curves show the two cuts ^(r, 0) and g{0,r), respectively. 

the figures at large x, while the total energies per particle are almost the same when prop- 
erly scaled with the scattering length. This points towards a delicate balance between the 
kinetic and potential contributions, which change with a but keep their sum constant once 
expressed in scattering length units. 

The last quantity analyzed in this work is the one-body density matrix pi{ri,r'j), which 
provides a measure of the overlap between two instances of the ground state wave function 
when one particle is shifted from its initial position at ri to a new position at r'^ 

In the case of translationally invariant systems as the one under study, the one-body density 
matrix depends on its arguments only through their difference and thus pi(ri, r'^^) = pi(ri — 
r[, 0) = pi(rii/) Additionally, if the interaction is isotropic, pi depends only on the magnitude 
of its argument rn/ = |rii/| and its large-rn/ limit measures directly the condensate fraction 
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FIG. 7: (Color online) Static structure function for polarization angles a = 0.4 and a = 0.6 
for X = 10"^, 10"'^ and 10~^. The red and blue curves show the two cuts S{k,0) and S{0,k), 
respectively. 

no which is proportional to the number of particles in the Bose-Einstein condensate. In the 
present case, however, the system is homogeneous but not isotropic so pi(rii/) will depend 
on the direction of rn/. Due to translational invariance, though, momentum is still a good 
quantum number and one expects condensation to appear at the zero momentum state. In 
that sense one can still write the relation between pi(rii/) and the momentum distribution 
in the form 

Pi(riiO = pno + J rfke^^"'n(k) (27) 

where n(k) is the momentum distribution of the non-condensate atoms. The one-body 
density matrix of the anisotropic gas of Bose dipoles can be further expanded in partial 
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FIG. 8: (Color online) Cuts of the one-body density matrix along the x (red dotted lines) and y 
(blue lines) axes, for the gas parameter values x = 10~^ and x = 10^^ (top to bottom). The curves 
on the left and right panels correspond to a = 0.4 and a = 0.6, respectively. 



waves 

oo 

Pi (r) = pim (r) cos(2m 6) , (28) 

m=0 

with pim{f) the radial function corresponding to the m-th mode contribution. Notice that, 
as before, the Bose symmetry restricts the previous sum to even modes only. 

Once enough modes pim{f) are known, one can reconstruct the complete one-body density 
matrix for all points in the plane. In particular, the cuts along the two directions parallel 
and perpendicular to the polarization plane, corresponding to = and 7r/2 in Eq. fl28l) . 
turn out to be particularly easy to evaluate 

oo oo 

p,{r,0) = J2pUr) , Pi{0,r) = J2i-^rpim{r) (29) 

m=0 m=0 

and display the maximum difference two cuts along different directions can take at the low 
densities considered in this work. Figure [S] shows the parallel and perpendicular cuts of pi(r) 
for the polarization angles a = 0.4 and a = 0.6 (left and right panels). The upper and lower 
curves correspond to the gas parameter values x = 10~^ and x = 10~\ respectively. As 
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before, the coordinates on the horizontal axis have been scaled with the density. Similarly 
to what happens to the other quantities analyzed, only at the highest gas parameter values 
the effects of the anisotropy start to be visible. This stresses once again the minor role 
played by the anisotropy at low densities, even in a non-diagonal quantity like pi(rii/). 

The most significant differences in the one-body density matrix for different values of the 
gas parameter appear at large distances, where pi(rii/) reaches an asymptotic value that 
can be identified with the condensate fraction uq in isotropic systems. When the anisotropic 
character of the interaction is taken into account, the presence of higher order partial waves 
in Eqs. (1251) and (^U^ could in principle change this behavior, making the limiting value of 
Pi(rii/) depend on the direction. The role of the different partial waves in that limit can be 
determined by looking at the momentum distribution of the system, which can be obtained 
from pi(rii/) by looking at the inverse of Eq. ([2] 



'•oo /'27r 

n(k) = / drr ^ikrco^ie-^) 



JO 



+ ^Pim(^)cos(2m6') 



m=l 



(30) 



with (f the angle formed by k and the x-axis. Changing variables a = 6 — (f, using the 
Jacobi-Anger expansion of a plane wave in Bessel functions 

oo 

^ikzcosa ^ j^^f^^^ + 2 ^ i'^Jmikz) cos(ma) (31) 

m=l 

and taking into account the orthogonality of the cosine functions in the range [0,27r], one 
finally finds 

/•oo °° /"OO 

fi(k) = 27r / Jo(/i;r)[pio(r) - nojrcir + 27r V'(-l)™cos(2mv?) / J2m{kr)pi,n{r) rdr 

Jo ^ ^ m=l "^0 

(32) 

where the first term on the right is isotropic and constitutes the m = mode of n(k), while 
the other terms stand for the m > contributions. Notice once again that only even modes 
appear in this expansion. 

Requiring n(k) to be finite for all values of k implies all integrals appearing in Eq. fl32|) 
to be finite, a constraint that can only be fulfilled when the functions multiplying the Bessel 
functions decay to zero at large distances. This condition particularly means that no can be 
obtained as the large r limit of the m = mode of the one-body density matrix, which is the 
isotropic contribution to pi(rii/). This is the direct generalization of the usual procedure 
employed to determine no in homogeneous and isotropic systems. 
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FIG. 9: (Color online) Left panel: isotropic (m = 0) contribution to the one-body density matrix 
at = 0.1 for the three polarization angles a = 0.2,0.4 and 0.6 (red solid, blue dashed and black 
dotted lines, respectively). Right panel: Condensate fraction no as a function of the gas parameter 
for a = 0.2,0.4 and 0.6, compared with the Bogoliubov prediction (black line). The color coding 
for the symbols is the same as in the left panel. 

Figure [9] shows on the left panel the m = mode contribution pio{r) for the three 
polarization angles a = 0.2,0.4 and 0.6 in terms of the scaled distances n^/^r for x = 0.1. 
As can be seen from the figure, all three curves are hardly distinguishable, stressing once 
again that to a large extent the physics is governed by the scattering length, which makes 
the density change for different polarization angles when x is fixed. The right panel in the 
figure shows the condensate fraction as a function of the gas parameter x = no?, obtained 
from the r — )■ oo limit of a fit to the long range asymptotic limit of the m = partial wave 
contribution of the one-body density matrix. Up to the highest value of x considered all 
three cases yield nearly the same prediction within statistical errors, while differences start 
to be significant only at x ~ 0.1. Therefore, the scaling on the gas parameter is preserved 
although moving from a; = 0.2toQ; = 0.6 for fixed x implies a change in density by almost 
a factor of 2. The figure also shows the Bogoliubov prediction for an isotropic gas of weakly 
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interacting 2D bosons 

= 1 - ^ (33) 

which agrees reasonably well with the Monte Carlo prediction up to a; ~ 0.01 where particle 
correlations seem to deplete the condensate less effectively than the mean field model. 



IV. SUMMARY AND CONCLUSIONS 



To summarize, in this work we have described the ground state properties of a gas of fully 
polarized Bose dipoles moving on the XY plane, where the polarization field forms an angle 
a with the normal direction. The projection of the polarization vector on the XY plane 
defines the x-axis, where the potential is softer than in any other direction. In this context, 
the dipole-dipole interaction defines a critical angle ccc ~ 0.615 where the potential starts to 
have attractive contributions. We have solved the zero energy two-body scattering problem 
by means of a Green's function and a decomposition of the wave function in partial waves. 
We have then found the dependence of the s-wave scattering length on the polarization angle 
by inspection of the m = mode, which dominates at large distances. Equipped with the 
two-body solution, we have built a variational wave function of the Jastrow type that has 
been used as a guiding function in a DMC simulation of the gas of polarized dipoles at low 
densities. We have found that the scaling of the energy in the gas parameter is preserved up 
to values of x where other isotropic systems deviate significantly. This behavior extends to 
other relevant ground state quantities like the pair distribution function, the static structure 
factor and the one-body density matrix, including the condensate fraction which can be 
determined from the large distance asymptotic behavior of its isotropic part. 
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